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基于 蒙特 卡 罗 仿 真 的 湖 库 水 质 预测 及 富 营养 化 风险 评估 方法 
EDE, AFF, BRR, E V WEP, FAR 


(北京 工商 大 学 计算 机 与 信息 工程 学 院 , 北京 100048) 


摘 E: 湖 库 水 质 预 测 及 富 营养 化 分 析 是 水 污染 防 控 的 重要 技术 手段 。 然 而 ， 已 有 的 水 质 预测 研究 通常 是 单 值 预测 ， 
并 以 此 为 依据 分 析 富 营养 化 状态 ， 具 有 一 定 的 偶然 性 和 不 确定 性 。 结 合 水 质 动力 学 模型 ， 提 出 了 一 种 基于 蒙特 卡 罗 
仿真 的 湖 库 水 质 预 测 及 富 营养 化 风险 评估 方法 。 在 已 知 水 质 动力 学 模型 水 质 指标 和 模型 参数 的 先 验 分 布 基础 上 ， 利 
用 蒙特 卡 罗 仿 真 预测 水 质 指 标的 演化 过 程 ， 获 得 未 来 时 刻 水 质 指 标 取 值 的 概率 分 布 ， 实 现 水 质 预测 。 进 一 步 ， 构 造 
综合 营养 状态 指数 ， 结 合 水 质 指 标 预测 结果 ， 计 算 综合 营养 状态 指数 的 概率 分 布 和 处 于 不 同 营养 程度 的 概率 ， 实 现 
富 营养 化 风险 评估 。 仿 真 结 果 表明 ， 该 方法 能 够 有 效 实 现 水 质 预 测 和 富 营养 化 分 析 ， 且 考虑 更 加 全 面 、 准 确 ， 克 服 
了 单 值 预测 结果 带 来 的 偶然 性 。 
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Monte-carlo simulation based water quality prediction and 
eutrophication risk assessment method of lakes 


Wang Xiaoyi, Zhou Yuqin, Zhao Zhiyao, Wang Li, Xu Jiping, Yu Jiabin 
(School of Computer & Information Engineering, Beijing Technology & Business University, Beijing 100048, China) 


Abstract: Water quality prediction and eutrophication analysis of lakes are important technical means of water pollution 


prevention and control. However, the existing water quality prediction research is usually in a form of single-valued 
prediction, and analyze eutrophication status on this basis, which has a certain degree of haphazard and uncertainty. 
Combining with the water quality kinetic model, proposed a water quality prediction and eutrophication risk assessment 
method based on Monte-Carlo simulation. Based on the prior distribution of water quality index and model parameters of 
water quality kinetic model, used Monte Carlo simulation to predict the evolution of water quality index to obtain the 


probability distribution of water quality indicators in future time and achieve water quality prediction. Further, constructed 


an integrated eutrophication status index. Combining with the predicted results of water quality indexes, calculated the 
mu probability distribution of comprehensive nutritional status index and the probability of different nutritional status to assess 
the eutrophication risk. The simulation results show that the proposed method can effectively predict the water quality and 
analyze eutrophication status, with more comprehensive consideration and accuracy. Meanwhile, it overcomes the 
haphazard brought by single-valued prediction result. 
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0 引言 太湖 以 及 梨 湖 都 已 经 完全 变 为 富 营 养 化 状态 UH。 因 此, 水 环 
i 境 问 题 已 迫 不 容 缓 。 水 质 预 测 能 有 效 的 掌握 水 质 情 况 ， 预 防 

随 着 现代 社会 的 飞速 发 展 ， 水 体 富 营养 化 已 成 为 一 个 全 ”大 规模 水 污染 事件 的 发 生 ， 对 水 污染 防治 具有 现实 意义 。 所 

球 性 的 重大 水 环境 问题 。 目 前 ， 在 世界 范围 内 ， 受 富 营养 化 ”以 实现 合理 有 效 的 水 质 预测 ， 并 在 此 基础 上 评估 富 营 养 化 风 
影响 的 沿海 地 区 有 500 多 个 。 此 外 ， 许 多 世界 淡水 湖泊 、 溪 险 是 解决 水 环境 问题 的 重要 思路 。 


流 和 水 库 正 在 遭受 富 营 养 化 。 相 关 研 究 表 明 ， 我 国 当前 已 有 己 有 的 水 质 预测 方法 主要 分 为 两 类 : 
80% 的 湖泊 出 现 富 营 养 化 现象 。 在 我 国 除 一 些 人 烟 稀 少 的 淡 a) 基 于 机 理 的 水 质 预 测 方法 。 这 类 方法 主要 是 通过 分 析 
水 湖区 域 以 及 还 保持 着 原始 状态 的 湖泊 之 外 ， 其 他 湖泊 营养 。 ”水 体内 部 的 机 理 演化 过 程 来 建立 水 质 模型 ,再 基于 该 水 质 模 


盐水 平 的 浓度 已 经 超过 了 湖泊 富 营 养 化 的 标准 浓度 。 尤 其 是 。 型 对 水 质 指 标 进 行 预测 。 常 用 的 机 理 模 型 包括 : QUAL2E、 
我 国 的 几 大 淡水 湖 的 营养 盐 浓度 都 已 经 超过 了 磷 、 氮 富 营养 WASP 和 EFDC 等 。QUAL2E 是 一 个 一 维 水 质 模型 ， 适 用 于 
化 的 发 生 浓 度 ,， 总 氮 浓 度 甚 至 超过 10 倍 之 上 。 当 下 ,我国 的 模拟 完全 混合 的 枝 状 河流 水 质 ， 它 假定 存在 主流 输送 机 理 ， 
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王 小 艺 


即 假定 平流 与 扩散 混合 都 是 沿 着 河流 的 主流 向 ， 而 在 河流 的 


横向 与 垂 向 上 水 质 组 分 是 完全 均匀 混合 的 ， 允 许 河流 治 程 有 
多 个 污染 源 、 取 水 口 以 及 文 流 汇 入 ;其 基本 方程 是 一 维 平 流 


sim 


用 于 


E 


M. 


染 物 和 


平移 一 志 


广 散 质量 迁移 方程 5-9 。 
ulation program) 是 美国 环保 
各 类 
河流 一 维 


WASP(water quality analysis 
局 (EPA) 开 发 并 推荐 使 用 的 
水 体 水 质 模拟 与 分 析 的 软件 ， 它 可 用 于 模拟 水 动力 
不 稳定 流 、 湖 泊 和 河口 三 维 不 稳定 流 、 常 规 污 
9 毒 污染 物 在 水 中 的 迁移 和 转换 规律 ， 其 基本 原理 是 
扩散 方程 & 31。EFDC (environmental fluid dynamics 


code) 模 型 
现 河流 、 
力学 、 水 质 模拟 、 

有 限 差分 模型 9。 
b) 数 据 驱动 的 水 质 预测 方法 。 


型 和 工 


是 涵盖 一 维 到 三 维 地 表 水 水 质 的 数学 模型， En sc 
水 库 、 湿 地 系统 、 河 口 和 海洋 等 水 体 的 水 动 
污染 oo 是 一 个 多 参数 


湖泊 、 


这 类 方法 利用 数据 分 析 模 
对 水 质 指标 的 实测 数据 进行 建 模 ， 再 通过 此 模型 对 


水 质 指标 进行 


络 、 时 间 序 列 模型 和 支持 向 量 相 


预测 。 常 用 的 数据 驱动 模型 包括 : 人 工 神经 网 
等 ， 该 类 模型 均 是 利用 不 同 


时 刻 的 水 质 指标 观测 数据 建立 表征 水 质 演化 过 程 的 数据 驱动 


模型 ， 


有 一 定 优势 ， 然 而 已 有 方法 大 者 


化 预测 
和 可 信 


在 偶然 性 和 不 确定 性 。 


在 此 基础 上 实现 水 质 过 程 预测 中 1。 两 类 方法 各 
得 到 的 是 水 质 指标 的 取 值 变 


c 


〈 单 值 预测 )， 其 结果 通常 具有 一 定 的 偶然 性 ， 精 确 度 
度 难以 令 人 信服 ， 在 此 基础 上 的 富 营养 化 评价 也 会 存 


本 文 为 了 弥补 单 值 预测 方法 带 来 的 不 足 ， 提 出 了 一 种 基 
于 蒙特 卡 罗 仿 真 的 水 质 预测 及 富 营养 化 风险 评估 方法 。 蒙 特 
FE 罗 仿 真是 指 利用 大 量 重复 随机 实验 来 获取 问题 数值 解 的 一 
种 仿真 。 本 文 基于 湖 库 水 质 指标 演化 机 理 ， 利 用 非 线 性 连续 
微分 方程 建立 湖 库 水 质 动力 学 模型 ， 进 而 给 定 水 质 指 标 和 模 


型 


参数 初 值 ， 利用 蒙特 卡 罗 方法 对 水 质 指标 注 
JU. 获得 未 来 时 刻 水 质 指标 取 值 的 概率 分 布 , 实现 水 质 预测 。 
进一步 ， 结 合 水 质 指标 预测 结果 ， 利 | 


化 过 程 进行 预 


综合 营养 状态 指数 法 


进行 评估 ， 获 得 综合 营养 状态 指数 的 概率 分 布 和 处 于 不 同 营 


养 程度 的 概率 ， 实 现 富 
水 质 
水 体 生态 系统 主要 包含 以 下 四 个 演化 过 程 : 


1 


过 程 、 


pe 


E: 


积 物 需 氧 量 。 
VERLA RE 


it DO 和 BOD 含量 
硝化 ， 反 硝化 ， 沉 积 ， 植 物 生长 ， 


浮游 植物 动力 学 过 程 、 磷 循环 和 氮 
质 指 标 描述 
(BOD)、 有 机 磷 (OP). 
氮 (NH-N) 和 硝酸 盐 氮 (NO3-N)C5 27, 


营养 化 风险 评估 。 


旨 标 演化 机 理 及 动力 学 模型 


溶解 氧 平衡 
循环 ， 可 由 八 个 水 
S: AEA (DO)、 浮 游 植 物 (Phyt)、 碳 生化 需 氧 
正 磷酸 盐 PO) ANA (ON). 
溶解 氧 平衡 过 程 通 
过 程 包括 复 氧 ， 碳 质 氧化 ， 
Phyt 死亡 ，Phyt 呼吸 和 沉 


变化 体现 ， 


浮游 植物 动力 学 过 程 的 主要 研究 对 象 是 水 体 中 
(浮游 植物 是 一 个 生态 学 概念 ， 是 指 在 水 中 


AE 
ri B 


营 浮游 生活 的 微小 植物 。 通 常 浮游 植物 就 是 指 浮游 藻 类 ), 通 


过 Phyt 含量 变 
亡 和 沉降 。 
过 OP 和 
正 磷酸 盐 ) 转换 为 浮游 植物 磷 ， 
可 机 磷 ， 矿 化 作用 ， 以 及 沉降 吸附 作用 。 
形态 相互 转 j 


2a N 


H 


$i 


物 氮 转 换 为 无 机 氮 和 
以 及 沉降 作用 。 


H 


zb 


量变 化 体现 ， 


化 体现 , 过 程 包含 三 个 部 分 : 浮游 植物 的 生长 、 
磷 循 环 是 磷 的 三 种 物质 形态 相互 转换 的 过 程 ， 

PO4 含 量变 化 体现 ， 过 程 包含 四 个 部 分 : 无 机 磷 
浮游 植物 磷 转 换 为 无 机 磷 
氮 循 环 是 氮 的 四 
和 响 的 过 程 ， 通 过 ON、NH3-N 和 NO3-N 含 
过 程 包含 : 无 机 氮 转 换 为 浮游 植物 氮 ， 浮 游 植 


机 氮 , 矿 化 作用 ,硝化 和 反 硝 化 作用 
水 体 生态 系统 演化 过 程 记 -5 如 图 1 所 示 。 
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矿 化 Y 
| B NO;-N - T i OS 
NH.N Suc i »| Phyt Y >| C:N:P 
Y | 死亡 + 呼吸 | 
光合 T 吸收 
BOD po, et or 
呼吸 Y 
L DO = 7 
CO; f ~i 
图 1 水 体 生态 系统 演化 过 程 
Fig. 1 Water ecosystem evolution process 


基于 上 述 水 质 指标 演化 机 理 ， 水 质 动力 学 模型 如 式 (1) 


所 示 。 


{k 


其 中 : 
指标 DO、Phyt、BOD、OP、PO4、ON、NH3-N 和 NO3-N 的 
浓度 。 
1 ks, ks} 为 模 


^ C 
Cpo =K (C, — Cpo) -1.047K 20—)C, 
po = k(Cyo,. — Cpo) s "TM )Csop 


—4.9371k,( Coo )Cw_w 0.054k, 
2+Cpo 


+[(6.6667 + 4k3)ks —0.3333k,. ]C pry 
Cp, = [ks — ke —0. Un 


Crop = -[1.047k, i 


>+ 0.025k, ]Csop 


+ 2.6667k, Cp ME 
i 1Co 


Co, =0.0125(0.125 + ke)C pn 一 [0. MM Cows )+ ta- kio)]Cop 
lt+Cp 20 


Co, 7 [0.0125(0.125 + ke)—0.025ks]Cp +0. 2766 — (D 


Cor 


Co, =0.125(0.125 + Ke) Cp -1.08k( E Cw " 


k 
7390-7 50Co 

Cw, y 7 [0.125(0.125 + ke) —0.025kskis]C pry: 

E 


C prye 


+1.08ko i ín +1.08ks G ALD N 


Cus 108, 
24C 


)Cyg, -n E 0.25k; (L— ki3)C pry 
DO 


-1.08k;( )Cvwo_w 


c 


{Coo， Cii» Cron» Cop» Cro,» Con» CNHs_N» CNON } 分 别 Ed ZR 水 质 


Phyt 的 浓度 以 叶绿素 a (Chl_a) 的 浓度 为 准 ， 参 数 
型 参数 ， 其 含义 和 单位 如 表 1 所 示 。 


= T 
4 x- (Cpo; [07 了 Cgop , Cop; Cro, 了 Con H [o7 , Cuyo, ) 


且 假 设 模型 参数 是 慢 


0 — (kk, s, k, ks, Ks ke, Ks Ks kioskir» Ks kiz)" , 
变 的 ， 则 式 (1) 可 表示 为 
无 | [f(x,0) 
Pur © 
其 中 ， f(x,0) 是 水 质 动力 学 模型 函数 。 
需要 注意 的 是 ， 不 同 的 湖泊 、 河 口 等 水 体内 部 环境 有 差 
别 ， 水 体 特 征 也 不 一 样 ， 故 不 同 水 体 的 模型 参数 也 不 一 样 。 
为 使 水 质 模型 更 为 精确 、 人 合理， 针对 不 同 水 体 环境 时 ， 模 型 
参数 9 的 取 值 应 结合 实测 数据 〈 水 质 指 标的 实测 数据 ) 和 相 
关 参 数 率 定 方法 G139 确 定 。 
2 基于 蒙特 卡 罗 仿 真 的 水 质 预测 及 富 营养 化 风险 
评估 
2.1 蒙特 卡 罗 仿 真 
蒙特 卡 罗 仿 真 (Monte Carlo simulation) 是 指 利用 大 量 重 
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复 随 机 实验 来 获取 问题 数值 解 的 一 种 仿真 方法 。 基 本 原理 是 
当 问 题 或 对 象 本 身 具 有 概率 特征 时 ， 可 以 用 计算 机 模拟 的 方 
法 产生 抽样 结果 ， 根 据 抽样 计算 统计 量 或 者 参数 的 值 ， 随 着 
模拟 次 数 的 增多 ， 可 以 通过 对 各 次 统计 量 或 参数 的 估计 值 求 
平均 的 方法 得 到 稳定 结论 。 该 方法 基于 大 数 定律 (事件 A 发 
生 的 概率 可 以 近似 为 在 大 量 实验 中 事件 A 发 生 的 频率 )， 通 
常用 于 解决 优化 、 积 分 、 随 机 概率 分 布 等 数学 问题 中 解析 解 
无 法 精确 给 出 的 情景 扩 。 蒙 特 卡 罗 仿 真 方法 的 步骤 如 表 2 


H 


Al 水 质 动力 学 模型 参数 
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Table 1 Water quality mechanism model's parameters 
符号 含义 单位 
k 20'C 复 氧 速率 day” 
k 20C 还 原 速率 day" 
k, 20C 硝 化 速率 day“ 
ka 底 泥 需 氧 量 g/(m - day) 
k; 浮游 植物 生长 速率 day" 
ks 浮游 植物 死亡 速率 day 
k 20C 反 硝化 速率 day“ 
ks 有 机 物 沉 降 速率 míday 
ko 溶解 有 机 氮 矿 化 速率 day" 
kio 有 机 磷 溶 解 比例 None 
ki 有 机 毛 溶 解 比例 None 
k 浮游 植物 呼吸 温度 系数 None 
ks 所 优先 选择 系数 None 
表 2 蒙特 卡 罗 仿 真 方法 的 步 又 


Table 2 Steps of the Monte Carlo simulation method 


步 又 内 容 


构造 随机 模型 
1 D 根据 问题 和 实测 数据 构造 随机 模型 ; 
2) 问题 的 解 对 应 该 模型 中 随机 变量 的 某 些 特征 。 
产生 随机 数 
1) 根据 模型 中 各 个 随机 变量 的 分 布 产 生 随机 数 。 
随机 抽样 
3 D 根据 随机 变量 的 分 布 特性 ， 设 计 和 选取 合适 的 抽样 方法 ; 
2) 对 每 个 随机 变量 进行 随机 抽样 。 
仿真 实验 
1) 按照 所 建立 模型 进行 仿真 实验 、 计算, 求 出 问题 的 随机 解 。 
统计 分 析 
5 D 统计 分 析 模 拟 实验 结果 ， 给 出 问题 的 概率 解 以 及 解 的 精度 
估计 。 


22 基于 蒙特 卡 罗 仿 真 的 水 质 指标 演化 预测 
考虑 到 蒙特 卡 罗 仿 真 的 优势 ， 本 文 基于 蒙特 卡 罗 仿 真 对 
水 质 指标 演化 进行 预测 。 
首先 ， 式 (2) 可 写成 离散 化 形式 如 下 : 
x, 7x, HAX f (x, 4,0, 1) - 0,4, 
0, 20, * hx f (x, 4,0.) * Wo 


3) 


8xl 13x1 8xl 13x1 8xl 
号 中: x R*.9 R"U,f:R"xR" >R”, 


水 质 指标 
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的 过 程 噪声 w —.(0,2,) ， 模 型 参数 过 程 噪声 ws ~ A/(0.2,) ， 
其 中 二 、 杞 为 相应 的 噪声 协 方差 阵 ，hh<1day) 是 预测 步 长 。 


在 基于 蒙特 卡 罗 仿 真 的 预测 过 程 中 ， 水 质 指标 演化 过 程 
会 被 多 次 模拟 。 对 每 一 次 模拟 ， 每 一 个 水 质 指 标 在 每 一 个 预 
测 时 间 点 的 数值 均 不 相同 ， 这 些 数值 同时 具有 规律 性 和 随机 
性 。 规 律 性 体现 在 预测 值 是 基于 确定 的 水 质 动力 学 模型 、 水 
质 指标 和 模型 参数 初 值 而 产生 的 ， 随 机 性 体现 在 每 一 个 水 质 
指标 在 每 一 个 预测 时 间 点 受 噪声 影响 而 得 到 的 预测 值 不 尽 相 
同 。 当 模拟 次 数 充 分 大 时 ， 结 合 概率 统计 规律 ， 可 以 获得 
水 质 指标 在 不 同时 间 点 上 取 值 的 概率 分 布 。 具 体 流程 步骤 为 ， 


a) 令 水 质 指标 的 初 值 为 Xx。， 模 型 参数 的 初 值 为 ,蒙特 


卡 罗 仿 真 粒子 数 为 NN， 预测 天 数 为 T7， 预 测 步 长 为 
h X(hsldap , HJR ETZI] k=0 。 


b) 根 据 水 质 指标 和 模型 参数 的 初 值 Xo) op 个 粒 
FRYE, BI Vie[h,N] ， 令 

Xo = xX0,0" —6, (4) 

c) 对 于 时 刻 k=h,24,…,7 ,基于 k-1 时 刻 的 水 质 指标 和 模型 

参数 值 09.9.) ， 利 用 式 (3) 进 行 单 步 预测 ， 得 到 时刻 的 

水 质 指标 和 模型 参数 预测 值 00,6, ) 。 

d) 若 k<T，k < kk+h， 转 至 c)， 否 则 转 至 e). 


e) 对 vke[h,T] ， 计 算 水 质 指标 和 模型 参数 的 概率 密度 函 
数 ; 


PG, 6.8) YANS, -xD) 


a (5) 
P(8,65.6,) ~ > 0/ N)ó(8, - 8?) 


其 中 : Pale) 为 基于 (x9. 05) 的 水 质 指标 概率 密度 函数 ， 


POED) 为 基于 (xb) 的 模型 参数 概率 密度 函数 , 5 为 狄 拉 
2.3” 富 营养 化 风险 评估 
本 节 的 富 营养 化 风险 评估 方法 ， 基 于 上 述 水 质 指标 演化 


结果 ， 构 造 综合 营养 状态 指数 7Z7(Z) ,计算 70G 的 概率 分 


布 , 通过 对 TLS) 的 分 级 处 理 计算 处 于 不 同 营养 程度 的 概率 。 

选取 Chl_a〈 即 本 文 的 Phyt)、 总 磷 (TP) 和 总 氮 (TN) 三 个 
水 质 指标 ， 计 算 三 个 水 质 指标 的 营养 状态 指数 ， 通 过 权 值 归 
一 化 得 出 TLE); X} TUS) 做 分 级 处 理 ， 并 计算 得 出 TUE) 的 
概率 分 布 和 处 于 不 同 营 养 程度 的 概率 8-431。 有 具体 流程 步 又 
为 : 


a) 对 于 时 刻 上 = 及 24,…,T ， 粒 子 x6. Pell N] ; 


(a) 对 于 Chl a，Crwx eX? ， 有 
TLI, (Chl _ a)-10Q.500«-1.086In(C,,, , )) (6) 
XIF TP, Con Cro €x, 
TLI, (TP)=10[9.436+1.6241n(Cop x ES Cro, )] (T) 
对 于 TN, Cono Cyu nas Crone € Xi, 


TLI, (TN)=10[5.453+1.6941n(Coy 
+ Cyg, vx + Cuo, v 4)] 


(b) EA Chl, a 为 基准 参数 ， 则 Chl a. TP. TN 的 权重 分 别 


3 
Wj UNE (9) 
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HB. om. om 2326 Chl. a 的 权重 和 与 自身 的 相关 系数 ，w、 


72 TALE TP. 的 权重 和 与 基准 参数 Chl a 的 相关 系数 , ws、 
分 别 是 TN 的 权重 与 基准 参数 Chl. a 的 相关 系数 。 
结合 各 水 质 指 标的 营养 状态 指数 和 权重 ， 计 算 TLOD : 
TLI, (£) = w,-TLI, (Chl | a) 
w, *TLI,(TP) + w,-TLI,(TN) 
b) 对 TLUS) 做 分 级 处 理 。 TUSD 的 值 越 大 营养 程度 越 重 ， 
以 0 一 100 分 级 ， TLI) <30 为 贫 营养 (oligotrophic)， 


(10) 


pup 
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x, —[6.5,0.05019,3,0.07,0.06,0.35,0.36,0.1] , 水 质 指标 过 程 噪声 协 方 
差 阵 X, = diag(0.1,0.001,0.01,0.001,0.001,0.01,0.0,0.01) ， 模 型 参数 初 


， 模 型 参数 过 程 噪声 协 


富 营养 化 风险 评估 方法 


值 0 =[1.208,0.01,0.01,0.15,0.13,0.0009,0.001, 
0.025,0.075,0.9,0.35.1.02,1.3]" 


2, = diag(0.1,0.01,0.001,0.1,0.001,0.001,0.0001, 
在 “9 FE] 5H 3 
方差 阵 0.0001, 0.0001, 0.01,0.001,0.01,0.1} 。 中 国 湖 (水 


库 ) 的 基准 参数 Chl a 与 其 他 参数 之 间 的 相关 系数 i 如 下 表 3 
所 示 。 预 测 步 长 4&=1/24( 单 位 : day) 、 粒 子 数 N=1000 和 仿真 时 


30< TLI(Z) «50 为 中 营养 (mesotrophic), TLIG) >50 为 富 营 养 
(eutrophic). 
c) 计 算 TLIE) 的 概率 密度 函数 : 


PLL Œ) ~ > IQ/ N)óQLI, (3) -TL (E)9)] (11) 


dd) 根据 综合 营养 状态 指数 7TLC) 的 概率 分 布 计算 处 于 不 
同 营养 程度 的 概率 Pe, TUO) 。 


ps, (TLI, C2) "ES Ead (12) 

其 中 : Pe CLO) Jg gp Za] k N £i R H LLO 属于 区 间 

EH 的 概率 ， BLE,， Ena] 73 TLL GO) 属于 区 间 LE, E] 的 个 
"" [X [R] LE, ANNUM 的 分 级 区 间 。 


3 ”案例 验证 


给 定 水 质 指标 和 模型 参数 初 值 ， 基 于 MATLAB 验证 水 
质 预测 及 富 营养 化 风险 评估 方法 的 有 效 性 。 
3.1 仿真 配置 

本 文 仿真 配置 的 水 质 指标 初 值 


长 T=100( 单 位 : day, 
表 3 Chla 5 TP, TN 的 相关 关系 
Table 3 Correlation coefficients between Chl a and TP, TN 


Chl a TP TN 
n; 1 0.84 0.82 
rà 1 0.7056 0.6724 
w; 0.4205 0.2967 0.2828 
3.2 结果 与 分 析 


基于 上 一 小 节 的 仿真 配置 ， 利 用 3.2 节 中 给 出 的 水 质 指 
标 演 化 预测 算法 , 可 得 出 N 个 粒子 的 演化 曲线 , 如 图 2 所 示 ， 
其 中 和 轴 表 示 预 测 时 间 《〈 预 测 天 数 100 RK, 一 天 取 24 小 时 ， 
则 预测 时 间 为 2400 小 时 )，Y 轴 表 示 水 质 指标 的 取 值 范围 。 
为 进一步 清晰 地 表示 结果 ， 粒 子 形成 的 概率 分 布 的 概率 密度 
函数 如 图 3 所 示 ， 其 中 X 轴 为 预测 时 间 (为 使 结果 更 清晰 ， 


选取 预测 天 数 100 天 中 每 隔 10 天 的 同一 时 刻 为 预测 时 间 )， 
Y 轴 为 水 质 指标 的 取 值 范 围 ，Z 轴 为 水 质 指标 的 概率 密度 函 
数 。 


N 


PO, ( 单位 : mgP/L) 


NHN (单位 : mgN/L) 


70 80 90 100 


30 40 50 60 
时 间 (单位 ， 天) 


ON (单位 : mgN/L) 


NO,.N (单位 : mgN/L) 


70 80 90 100 


30 40 50 60 
时 间 (单位 ， 天) 


图 2 基于 蒙特 卡 罗 模 拟 的 水 质 预测 


Fig.2 Water quality prediction based on Monte Carlo simulation 
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EC 1- 
M Ras- 
E ges 4 0- 

0.08 


时 间 ( 单 位: 天) 


E 100 - : 
位 : 天 ) 时 间 (单位 ， 天) 
OP (单位 : mgP/L) 
x 1 x 1 
Bis. Ps 
ES Deu E 0- 
92 = 04 B 
100 E rue --— “100 
c o 
时 间 (单位 :天 ) a Fi hr. 
PO, Cfi: mgPiLO 时 间 (单位 :天 ) 
g 
Sg 
Eas 4 
其 | 
$ ol 
0.7 =h 
100 
° TE CE: 天 ) i 时 间 (单位 :天 ) 
NH,-N (单位 : mgN/L) NO,-N (单位 : mgN/L) 


图 3 特定 时 间 水 质 指 标的 概率 分 布 函数 


Fig.3 Probability distribution functions of water quality indicators at specific time indexes 


概率 密度 函数 


综合 富 营养 化 指数 90 0 时 间 〈 单 位 ; 天 ) 


图 4 TLIE) 在 特定 时 间 的 概率 分 布 函数 
Fig.4 Probability distribution functions of TL/(Z) at specific time indexes 


1 


08r 


02r 


0o 10 20 30 40 30 60 70 80 90 100 
时 间 (单位 ， 天) 
图 5 不 同 富 营 养 化 等 级 概率 曲线 


Fig.5 The probability curves of being in different eutrophication levels 
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区 
预测 值 的 取 值 范围 变化 大 , 发 散 怕 
其 中 
性 强 ， 水 质 指 标 
水 体 含 量 
围 内 。 这 也 体现 在 图 
随 着 预测 时 间 的 
EJE DRE IE AE 7 E Di 
并 受 噪声 影响 ， 


H PO4 的 


Esth, 水 质 预 测 


线 可 知 
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所 示 的 8 个 水 质 指标 预测 值 


， 水 质 指标 


有 随机 性 。 


先 呈 现 较 快 的 下 降 趋 势 ， 再 


OP 和 PO4 的 发 散 性 较 小 ， 水质 


指标 DO 的 
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风险 评估 方法 考虑 到 水 质 预 测 的 随机 性 ， 弥 补 了 单一 定 值 的 
水 质 预 测 和 富 营 养 化 评价 方法 的 不 精确 性 ， 降 低 了 预测 的 偶 
然 性 ， 提 高 了 水 体 的 富 营养 化 风险 评估 的 准确 度 。 为 进一步 


提高 预测 的 精度 ， 后 期 工作 主要 分 为 两 个 部 分 : 采用 参数 率 


稳定 在 较 


2 
取 

， 水 质 指 标 Phyt, BOD, ON, NHs-N 和 NO3-N 的 发 散 
首 
也 


3 的 水 质 指标 预测 值 概率 分 


cr 


小 的 取 值 范 
布 曲线 中 。 


曾 加 ，8 个 水 质 指标 的 预测 值 概 


率 分 布 由 访 


状 , 这 体现 了 水 质 指标 的 
对 噪声 敏感 。 其 中 ， 水 质 指标 Phyt、BOD、 
ON, NH-N 和 NOs-N 的 随机 性 较 强 ， 受 
声 敏 感 ， 水 质 指标 OP F 


标 预 测 值 
质 指 标 预测 


噪声 敏感 度 较 小 ; 
寸 噪声 敏感 度 小 。 
线 和 概率 分 布 


受 噪 声 影 响 ， 


2r 
综合 图 


具有 随机 性 ， 


Li —L E 
PA Jb 


2 和 3 所 示 的 


响 大 ， 对 噪 


随机 性 较 小 ， 受 噪声 影响 较 
水 质 指标 DO 的 随机 性 


小 ， 受 噪声 
8 个 水 质 指 


线 可 知 ， 水 体 环 境 较 为 复杂 ， 水 


范围 内 ， 上 有 具有 确定 性 。 无 儿 
变化 速度 较 慢 。 


界 特 殊 


AR 


向 时 ， 


有 随机 性 ， 但 预测 值 稳定 在 一 定 


此 阶段 水 体 


EFE 3 给 出 的 水 质 指标 概率 分 布 ， 利 用 3.3 节 给 出 的 
富 营 养 化 风险 评估 方法 ， 可 得 到 综合 营养 状态 指数 Tu OD 的 
Y- 概率 分 布 ， 以 及 处 于 不 同 营养 状态 的 概率 ， 结 果 如 图 4、5 
Z 。 ”所 示 。 图 4 中 ,，X 轴 表 示 预 测 时 间 (为 使 结果 更 清晰 ， 选 取 
Cip 预测 天 数 100 天 中 每 隔 10 天 的 同一 时 刻 为 预测 时 间 )，Y A 
= 表示 综合 营养 状态 指数 的 取 值 范围 ，Z 轴 为 综合 营养 状态 指 
O 数 的 概率 密度 函数 。 图 5 中 ，X 轴 表 示 预 测 时 间 (预测 天 数 
© 100 天 ， 一 天 取 24 小 时 ， 则 预测 时 间 为 2400 小 时 )，Y 轴 表 
v0 示 不 同 营养 状态 的 概率 。 
e 图 4 中 ， 由 水 体 富 营 养 化 概率 密度 函数 可 知 富 营养 化 风 
C» 。 险 评估 具有 随机 性 。 随 着 预测 时 间 的 增加 ， 富 营养 化 风险 评 
e cote 4 PUETE FR AK, B e PECORE DIS CIAR 
Y 有 随机 性 和 偶然 性 。 这 也 体现 在 图 5 中 ， 随 着 预测 时 间 的 
>= ”增加 ， 贫 营养 状态 概率 减 小 ， 中 营养 状态 概率 增加 。 据 浮游 
会 ”植物 特性 可 知 ， 浮 游 植物 的 生长 具有 季节 性 ， 水 体 中 的 富 营 
Sé — 养 程度 一 般 也 具有 季节 性 。 本 文选 取 的 水 质 指标 和 模型 参数 
CD) 初 值 为 水 体 贫 营 养 状 态 ， 且 水 体 营养 物质 含量 为 增长 阶段 。 
C. 综 上 所 述 ， 基 于 水 质 预 测 的 水 体 富 营养 化 风险 评估 具有 规律 
co 性 和 随机 性 。 
O 综 上 所 述 ， 无 论 是 水 质 指标 预测 值 的 概率 分 布 ， 还 是 富 
营养 化 风险 评估 的 概率 分 布 可 知 ， 水 体 环境 由 多 因素 影响 ， 
了 复杂 性 ， 所 以 水 质 预测 具有 随机 性 。 本 文 提 出 的 基于 水 
质 动力 学 模型 和 蒙特 卡 罗 仿 真 的 水 质 机 理 预测 方法 是 概率 性 
预测 ， 相 较 于 单一 取 值 的 水 质 预测 ， 精 确 度 更 高 。 
4 ”结束 语 
本 文 结合 水 质 动力 学 模型 ， 提 出 了 一 种 基于 蒙特 卡 罗 仿 


真 的 湖 库 水 质 预 测 


指标 演化 机 理 ， 


利用 非 线 性 


及 富 营 养 化 风险 评估 方法 。 基 于 湖 库 水 质 


连续 微分 方程 建立 湖 库 水 质 动力 


学 模型 。 在 已 知 水 质 动力 学 模型 水 质 指标 和 模型 


参数 的 先 验 


分 布 基础 上 ， 利 用 蒙特 卡 罗 仿 真 预 测 8 个 水 质 指标 的 演化 过 
程 ， 获 得 8 个 水 质 指标 取 值 的 概率 分 布 ， 实 现 水 质 预测 ， 从 
仿真 结果 可 知 水 质 指标 预测 具有 随机 性 ， 且 不 同 水 质 指标 对 
噪声 的 敏感 程度 不 一 。 进 一 步 ， 本 文通 过 构造 综合 营养 状态 


指数 ， 结 合 水 质 指标 预 涡 


率 分 布 ， 同 时 ， 


| 结果， 计算 综合 营养 状 


态 指数 的 概 


对 综合 营养 状态 指数 分 级 处 理 ， 获 得 处 于 不 


同 营养 程度 的 概率 ,实现 富 营养 化 风险 讨 


基于 水 质 预 测 的 富 营养 化 风险 评估 


本 文 所 提出 的 基 


也 具有 随机 性 


EF 估 ; 仿真 结果 可 知 ， 


。 综 上 所 述 ， 


于 蒙特 卡 罗 仿 真 的 湖 库 水 质 预测 


及 富 营养 化 


定 方法 ， 结 合 水 质 指标 的 实测 数据 ， 获 得 水 质 动力 学 模型 的 
模型 参数 的 估计 值 ， 采 用 一 定 的 方法 克服 模型 参数 的 慢 变 假 


设 。 
参考 文献 : 


[1] Meinikmann K, Hupfe M r, Lewandowski J. Phosphorus in groundwater 


discharge: a potential source for lake eutrophication [J]. Journal of 
Hydrology, 2015, 524 (5) :214-226. 
[2] Rast W, Thornton J. Trends in eutrophication research and control [J]. 
Hydrological Processes, 2015, 10 (2): 295-313. 
[3] Paerl H, Xu Hai, Hall, N S, et al. Nutrient limitation dynamics 
examined on a multi-annual scale in Lake Taihu, China: implications 
for controlling eutrophication and harmful algal blooms [J]. Journal of 
Freshwater Ecology, 2015, 30(1): 5-24. 
TE, URF, 等 . 我 国 湖泊 富 营养 化 问题 及 其 控制 措施 研 
ZK [I]. 绿色 科技 , 2015, 2015(6): 189-191. (Wu Yaojie, Yu Xia, 


[4] X5, 
Liu Dingfu, et al. Research progress on eutrophication of lakes in China 
and its control measures [J]. Journal of Green Science and Technology, 
2015, 2015(6): 189-191. ) 

[5] Woldegiorgis B T, Van Griensven A, Pereira F, et al. A new 
unconditionally stable and consistent quasi-analytical in-stream water 
quality solution scheme for CSTR-basedwater quality simulators [J]. 
Water Resources Research, 2017, 53 (6): 4668-4690. 

朱文 婷 ， 钱 新 ， 钱 玮 ， 等 . 基于 QUALOK 模型 的 尾 水 生态 净化 方案 
设计 [D]. 环境 化 学 , 2018, 37(3): 600-608. (Zhu Wenting, Qian Xin, 
Qian Wei, et al. Application of QUAL2K model for ecological 


[6] 


purification design of treated water [J]. Environmental Chemistry, 2018, 
37(3): 600-608. ) 

HI, 卢 诚 , (RUE, 等 . 滥 水 河流 域 水 质 时 空 变化 特征 及 其 污染 源 解 
析 J]. 环境 科学 学 报 , 2017, 37(8): 2829-2837. (Qiu Yu, Lu Cheng, 


[7] 


Xu Zhe, et al. Spatio-temporal variation characteristics and water 
pollution sources in the Huangshui River Basin [J]. Acta Scientiae 
Circumstantiae, 2017, 37(8): 2829-2837. ) 

[8] Cho Jaepil, Her Younggu, Bosch D. Sensitivity of Simulated 
conservation practice effectiveness to representation of field and 
in-stream processes in the little river watershed [J]. Environmental 
Modeling & Assessment, 2016, 22(2): 1159-173. 

[9] Kim T, Sheng Y P. Estimation of water quality model parameters [J]. 

Ksce Journal of Civil Engineering, 2010, 14(3): 421-437. 

[10] 朱文 博 , 王 洪 秀 , MMA, SB. 河道 曝 气 提升 河流 水 质 的 WASP 模型 研 
5L 加. 环境 科学 ，2015，36(4): 1326-1331. (Zhu Wenbo, Wang 
Hongxiu, Liu Cui, et al. Improvement of river water quality by aeration: 
WASP model 2015, 36(4): 
1326-1331.) 


[11] Wang Xiaoyi, Jia Jie, Su Tingli, etal. A fusion water quality soft-sensing 


study  [J].Environmental Science, 


method based on WASP model and its application in water 

eutrophication evaluation [J]. Journal of Chemistry, 2018, 2018: 1-14. 
[12] Huang Jiacong, Gao Junfeng, Hórmann G. 
Hydrodynamic-phytoplankton model for short-term forecasts of 
phytoplankton in Lake Taihu, China [J]. Limnologica-Ecology and 
Management of Inland Waters, 2012, 42(1): 7-18. 


[13] Akomeah E, Chun K P, Karl-Erich Lindenschmidt. Dynamic water 


录用 定稿 王 小 艺 ， 等 : 


quality modelling and uncertainty analysis of phytoplankton and 
nutrient cycles for the upper South Saskatchewan River[J]. 
Environmental Science and Pollution Research, 2015, 22(22): 
18239-18251. 

[14] ERAK, HKR, hii, 等, 基于 EFDC 模型 对 好 阳 湖 湖面 面积 的 
推算 方法 [n 南京 信息 工程 大 学 学 报 ，2016，8(1):56-63.(Xiong 
Jiaging, Lai Geying, Peng Xiaojuan, etal. Poyang Lake area estimation 
based on EFDC model [J]. Journal of Nanjing University of Information 
Science and Technology: Natural Science Edition, 2016, 8(1): 56-63. ) 

[15] Arifin R R, James S C, De Alwis D A, et al. Simulating the thermal 
behavior in Lake Ontario using EFDC [J]. Journal of Great Lakes 
Research, 2016, 42(3): 511-523. 

[16] Zhang Chunxue, You Xueyi. Application of EFDC model to grading the 
eutrophic state of reservoir: case study in Tianjin Erwangzhuang 
Reservoir, China [J]. Engineering Applications of Computational Fluid 
Mechanics, 2017, 11(1): 111-126. 

[17] Gani Y, Ustinova E, Ajakan H, etal. Domain-adversarial training of 
neural networks [J]. Journal of Machine Learning Research, 2017, 17 
(1): 2096-2030. 

[18] 5L, B HH. 基于 灰色 模型 和 模糊 神经 网 络 的 综合 水 质 预测 模型 
研究 [J]. 环境 工程 学 报 ，2015，9(2): 537-545. (Zhang Ying, Gao 
Qianqian. Comprehensive prediction model of water quality based on 
grey model and fuzzy neural network[J].Chinese Journal of 
Environmental Engineering, 2015, 9(2): 537-545. ) 

[19] 5339, JM X, EAR. 基于 BP 神经 网 络 的 大 伙房 水 库 水 质 综合 
评价 D]. 沈阳 农业 大 学 学 报 ,，2014, 45(5): 637-640. Ma Lingling, 
Zhou Linfei, Wang Tieliang. Comprehensive evaluation of water quality 
in Dahuofang reservoir based on BP neural network [J]. Journal of 
Shenyang Agricultural University, 2014, 45(5): 637-640. 

[20] Yan Qiao, Ma Cong. Application of integrated ARIMA and RBF 
network for groundwater level forecasting [J]. Environmental Earth 
Sciences, 2016, 75(5): 1-13. 

[21] Gebler D, Wiegleb G,  Szoszkiewicz K. Integrating river 
hydromorphology and water quality into ecological status modelling by 
artificial neural networks [J]. Water Research, 2018, 139: 395-405. 

[22] Liu Binsheng, Wang Ying, Hu Xueping. A chaotic prediction method to 
time series data [J]. Advanced Materials Research, 2010,113-116 
(1662-8985): 1367-1370. 

[23] Katimon A, Shahid S, Mohsenipour M. Modeling water quality and 
hydrological variables using ARIMA: a case study of Johor River, 
Malaysia [J]. Sustainable Water Resources Management, 2017, 2017 
(03): 1-8. 

[24] Faruk D Ö. A hybrid neural network and ARIMA model for water 
quality time series prediction [J]. Engineering Applications of Artificial 
Intelligence, 2010, 23 (4): 586-594. 

[5] KEF, & R44, RE. 基于 人 工 鱼 群 优化 支持 向 量 机 水 文 预报 系 
统 模型 [J]. 计算 机 应 用 研究 ，2010，27(8): 2902-2905. (Zhang 
Ya-Ping, Lu Xian-Hua, Jin Mei. Hydrological forecasting system model 
based on AFSVM [J]. Application Research of Computers, 2010, 27(8): 
2902-2905. ) 

[26] Defne Z, Spitz F J, De Paul V, et al. Toward a comprehensive 
water-quality modeling of barnegat bay: development of ROMS to 
WASP coupler [J]. Journal of Coastal Research, 2017, 78: 34-45. 

[27] Jia Haifeng, Zhang Yansong, Guo Yu. The development of a 
multi-species algal ecodynamic model for urban surface water systems 


and its application[J]. Ecological Modelling, 2017, 221(15): 1831-1838. 


基于 蒙特 卡 罗 仿 真 的 湖 库 水 质 预测 及 富 营养 化 风险 评估 方法 


[28] 港 青青 . WASP 水 质 模型 和 GIS 空间 分 析 在 水 污染 监测 中 的 集成 应 
用 D]. 科技 资讯 ,2015, 13(34): 12-14. (Zhan Qingin. WASP water 
quality model and GIS spatial analysis integrated application in water 
pollution monitoring [J]. Science & Technology Information, 2015, 
13(34): 12-14. ) 

[29] Jtr, P, JA. 300 年 来 壮阳 湖 营 养 盐 演化 重建 与 模拟 D] E 
态 学 报 ，2016，36(8): 2393-2402. (Liao Mengna，Yu Ge，Guo Ya. 
Reconstructions and simulations of nutrient evolution in Poyang Lake 
over the past 300 years [J]. Acta Ecologica Sinica, 2016, 36(8): 
2393-2402. ) 

[30] Gooddy D C, Lapworth D J, Bennett S A., etal. A multi-stable isotope 
framework to understand eutrophication in aquatic ecosystems. [J]. 
Water Research, 2016, 88: 623-633. 

[31] Li Ranran, Zou Zhihong. Water environmental capacity analysis of 
taihu lake and parameter estimation based on the integration of the 
inverse method and Bayesian modeling [J]. International Journal of 
Environmental Research and Public Health, 2015,12(10): 12212-12224. 

[32] Shi Yuning, Davis K J, Zhang Fuging, et al. Parameter estimation of a 
physically-based land surface hydrologic model using an ensemble 
Kalman filter: A multivariate real-data experiment [J]. Advances in 
Water Resources, 2015, 83: 421-427. 

[33] KAŽ, EZ, Wih, $. WASP 模型 参数 率 定 与 敏感 性 分 析 [I]. 
水 资源 与 水 工程 学 报 , 2009, 20(5): 28-30. (Zhang Yongxiang, Wang 
Lei, Yao Weitao, et al. Calbration and sensitive analysis on the 
parameters of the WASP model [J]. Journal of Water Resources &W 
ater Engineering, 2009, 20 (5): 28-30. ) 

[34] Ala-Luhtala J, Whiteley N, Heine K, et al. An introduction to twisted 
particle filters and parameter estimation in non-linear state-space 
models [J]. IEEE Trans on Signal Processing, 2015, 64(18): 4875-4890. 

[35] He Wei, Williard N, Osterman M, et al. Prognostics of lithiumion 
batteries based on Dempster-Shafer theory and the Bayesian Monte 
Carlo method [J]. Journal of Power Sources, 2011, 196(23): 
10314-10321. 

[36] 朱 伦 伦 ， 能 九龙 , WEH, 等 . 基于 蒙特 卡 罗 仿 真 的 校准 结论 风险 分 
Ar [] 电子 测量 与 仪器 学 报 , 2016, 30(7): 1045-1051. (Zhu Lunlun, 
Xiong Jiulong，Xie Jinzhe，et al. Calibration conclusion Tisk analysis 
based on Monte Carlo simulation [J]. Journal of Electronic 
Measurement and Instrument, 2016, 30(7): 1045-1051. ) 

[37] Fougere N, Altwegg K, Berthelier J J, et al. Three-dimensional direct 
simulation Monte-Carlo modeling of the coma of comet 
67P/Churyumov-Gerasimenko observed by the VIRTIS and ROSINA 
instruments on board Rosetta [J]. Astronomy and Astrophysics, 2016, 
588 (A&A): A134. 

[38] & S, ZER, MER, SR. 基于 蒙特 卡 洛 小 波 去 噪 的 股票 投资 组 
合 风 险 优化 研究 [] 计算 机 应 用 研究 ，2018，35(10): 
2947-2951,3028. (Li Junchang, Fan Chongjun, Yang Yunpeng, et al. 
Stock portfolio risk optimization based on filtering method of wavelet 
using Monte Carlo simulation [J]. Application Research of Computers, 
2018, 35 (10): 29047-2951,3028.) 

[39] Chen Bin, Siepmann J I. A novel Monte Carlo algorithm for simulating 
strongly associating fluids: applications to water, hydrogen fluoride, 
and acetic acid [J]. Journal of Physical Chemistry B, 2015, 104 (36): 
8725-8734. 

[40] Jorgensen W L, Jenson C. Temperature dependence of TIP3P, SPC, and 
TIP4P water from NPT Monte Carlo simulations: Seeking temperatures 


of maximum density [J]. Journal of Computational Chemistry, 1998, 


201901.00024v1 


chinaXiv 


[41] 


[42] 


19(10): 1179-1186. 

李 镇 镇 ， FER, EK, F. 基于 动态 聚 类 分 析 和 盲 数 理论 的 综合 营 
养 状态 指数 评价 模型 [T]. 环境 工程 学 报 , 2015, 9(4): 2021-2026. (Li 
Zhenzhen, Li Xiaodong, Li Feil, etal. Improved assessment model for 
comprehensive trophic state index based on dynamic cluster analysis 
and blind theory [J]. Chinese Journal of Environmental Engineering, 
2015, 9(4): 2021-2026. ) 

MRA, Se, EDE. 基于 综合 营养 状态 指数 法 的 石 佛 地 水 库 水 质 
富 营养 化 评价 [I]. 沈阳 农业 大 学 学 报 , 2016, 47(1): 119-123. (GUO 
Cheng-jiu, HONG Mei, YAN Bin. Eutrophication of Shifosi Reservoir 


ChinaXiv 合 作 期 刊 


录用 定稿 王 小 艺 ， 等 : 基于 蒙特 卡 罗 仿 真 的 湖 库 水 质 预测 及 富 营养 化 风险 评估 方法 第 37 卷 第 3 期 


Based on Comprehensive Nutrition State Index [J]. Journal of Shenyang 
Agricultural University, 2016, 47(1): 119-123. ) 

[43] 李 林 衡 ， 郑 飞 ， 何 春花 ， 等 . PRO E UU JS CHE E RUE ORE 
营养 化 程度 D 中 国 给 水 排水 , 2016, 32(13): 75-78. (Li Linheng, 
Zheng Fei, He Chunhua, et al. Evaluation of reservoir eutrophication in 
Yinzhou district based on comprehensive eutrophication state index 


method [J]. China Water and Wastewater, 2016, 32(13): 75-78. ) 


